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Abstract 

We discuss how the main features of high-energy 'soft' and 'semihard' pp collisions may 
be described in terms of parton cascades and multi-Pomeron exchange. The interaction 
between Pomerons produces an effective infrared cutoff, fc sa t, by the absorption of low 
kt partons. This provides the possibility of extending the parton approach, used for 
'hard' processes, to also describe high-energy soft and semihard interactions. We outline 
a model which incorporates these features. Finally, we discuss what the most recent LHC 
measurements in the soft domain imply for the model. 



1 A unified description? 

'Soft' and 'hard' high-energy pp interactions are described in different ways. The appropriate 
formalism for high-energy soft interactions is based on Reggeon Field Theory with a phe- 
nomenological (soft) Pomeron, whereas for hard interactions we use a QCD partonic approach, 
where the (QCD) Pomeron is associated with the BFKL vacuum singularity pQ. However, the 
two approaches appear to merge naturally into one another. That is, the partonic approach 
seems to extend smoothly into the soft domain. 

The BFKL equation describes the development of the gluon shower as the momentum 
fraction, x, of the proton carried by the gluon decreases. That is, the evolution parameter is 
ln(l/x), rather than the ln/c^ evolution of the DGLAP equation. Formally, to justify the use 

1 Presented at Linear Collider 2011: Understanding QCD at Linear Colliders in searching for old and new 
physics, 12-16 September 2011, ECT*, Trento, Italy 



of perturbative QCD, the BFKL equation should be written for gluons with sufficiently large 
k t . However, it turns out that, after accounting for NLL(l/x) corrections and performing an 
all-order resummation of the main higher-order contributions [2], the intercept of the BFKL 
Pomeron depends only weakly on the scale for reasonably small scales. The intercept is seen to 
be A = Qjp(O) — 1 ~ 0.35 over a large interval of smallish k t , Fig. ??. Thus the BFKL Pomeron 
is a natural object to continue from the 'hard' domain into the 'soft' region. 



BFKL stabilized 



A = ap(0) - 1 



LL1/x: A = a s 4ln2 




resummed 
BFKL 



► a, 



Intercept A = a P (0) -1 ~ 0.35 
A depends weakly on k t 
for low k t 



NLL1/x: A = A (1-6.5 ot s ) 

Figure 1: The behaviour found for the Pomeron intercept at leading and next-leading log(l/x) 
order, where a s = a s /3ir. When an all-order resummation of the main high-order contributions is 
included, A tends to a value of about 0.35 for reasonably large values of a s . 



Figure 2: The cascade structure of a gluon ladder. The BFKL or QCD Pomeron is the sum of 
ladder diagrams, each with a different number of rungs. 



The BFKL or QCD Pomeron may be viewed as a sum of ladders based on the exchange 
of two t-channel (Reggeized) gluons. Each ladder produces a gluon cascade which develops in 
ln(l/x) space, and which is not strongly ordered in k t , see Fig. [2| There are phenomenological 
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arguments (such as the small slope of the Pomeron trajector)|^J the success of the Additive 
Quark Model relation^ etc.) which indicate that the size of an individual Pomeron is relatively 
small as compared to the size of a proton or pion etc. Thus we may regard the cascade as a 
small-size 'hot-spot' inside the colliding protons. 

At LHC energies the interval of BFKL ln(l/ x) evolution is much larger than that for DGLAP 
lnfc 2 evolution. Moreover, the data already give hints that we need contributions not ordered in 
k t , a la BFKL, since typically DGLAP overestimates the observed (k t ) and underestimates the 
mean multiplicity [H |S] . Further, it is not enough to have only one Pomeron ladder exchanged; 
we need to include multi-Pomeron exchanges. 

Basically, the picture is as follows. In the perturbative domain we have a single bare 'hard' 
Pomeron exchanged with a trajectory a h p TG ~ 1.35 + cn' hare t, where a' bare ^ 0.05 GeV -2 . The 
transition to the soft region is accompanied by absorptive multi-Pomeron effects, such that an 
effective 'soft' Pomeron may be approximated by a linear trajectory a e p ~ 1.08 + 0.25t in the 
limited energy range up to Tevatron energies [B] . This smooth transition from hard to soft is well 
illustrated by Fig. [3j which shows the behaviour of the data for vector meson (V = p, oj, 4>, J/ip) 
production at HERA, 7*p — > V(M)p, as Q 2 + M 2 decreases from about 50 GeV 2 towards zero. 



2 Multi-Pomeron diagrams 

The eikonal model accounts for the multiple rescattering of the incoming fast particles. We 
haveS 

ImT = (l-e- n/2 ) = (0/2) - (0 2 /8) + ... (1) 

which displays the multi-Pomeron corrections to the bare Pomeron amplitude, 0/2, that tame 
the power growth of the cross section with energy. Simultaneously, these multi-Pomeron dia- 
grams also explain the growth of the central plateau jU [3] 

dN diVi-Po m 
dr] ' drj ' 

where dNi_ Pom / dr] is the plateau due to the exchange of one Pomeron, which is independent 
of collider energy. The growth is due to the increasing number, np, of Pomerons exchanged 
as energy increases. These (eikonal) multi-Pomeron contributions are included in the present 



2 Recall that a' P oc l/(kf) oc Rp om - 



3 The argument is best seen by analogy with nuclear physics. For light nuclei we have 'additive' cross sections, 
a = A\Aio nrll since the nuclei radii R 3> r nn . On the other hand for a heavy nucleus, where r nn ~ R, large 
Glauber corrections break the additive result. Similarly, the experimental success of the AQM indicates that 

^ qq ^° Rpora Rp- 

4 To allow for low-mass proton dissociation, the amplitude ([l]) is written in matrix form, Tj&, between (Good- 
Walker [7]) diffractive eigenstates. 
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Figure 3: The parameters of the Pomeron trajectory, otp(t) = ap(0) + a' P t, determined from the 
energy and t behaviour of high energy HERA data for vector meson production, 7*p — > V(M)p. 



Monte Carlos to some extent, as a Multiple Interaction (MI) option, but Pomeron- Pomeron 
interactions are not allowed for. 

Since the (small size) Pomeron cascades (hot spots) occur at different impact parameters, 
b, there is practically no interference between them. Moreover, at this 'eikonal' stage, the 
multi-Pomeron vertices, which account for the interaction between Pomerons, are not yet in- 
cluded in the formalism. These are interactions between partons within an individual hot spot 
(Pomeron). Formally, these are NNLO interactions, but their contribution is enhanced by the 
large multiplicity of partons within a high-energy cascade. In terms of Reggeon Field The- 
ory, the additional interactions are described by so-called enhanced multi-Pomeron diagrams, 
whose contributions are controlled by triple-Pomeron (and more complicated multi-Pomeron) 
couplings^} Recall that non-enhanced (eikonal) multi-Pomeron interactions are caused mainly 
by Pomerons occurring at different impact parameters, and well separated from each other in 
the 6-plane. On the other hand, the enhanced contributions mainly correspond to additional 

5 These diagrams are responsible for high-mass proton dissociation. 
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interactions (absorption) within an individual hot spot, but with the partons well separated in 
rapidity. 

The main effect of the enhanced contribution is the absorption of low k t partons. Note that 
the probability of these additional interactions is proportional to <7 a b s ~ l/^t> an<J their main 
qualitative effect is to induce a splitting of low k t partons into a pair of partons each with lower 
x, but larger k t . Effectively this produces a dynamical infrared cut-off, k sat , on k t , and partly 
restores a DGLAP-like k t - ordering within the cascade at larger k t . 

3 Schematic sketches of the model 

Qualitatively, the structure of soft interactions based on the 'BFKL' multi-Pomeron approach is 
as follows. The evolution produces a parton cascade which occupies a relatively small domain in 
6-space, as compared to the size of the proton. We have called this a hot spot. The multiplicity 
of partons grows as x~ A , while the /c t 's of the partons are not strongly ordered and depend 
weakly on Ins. Recall A = op(0) — 1. Allowing for the running of a s , the partons tend to drift 
to lower kt where the coupling is larger. This is shown schematically in Fig. 111(a). 



(a) single BFKL Pom. (b) DGLAP-based MC (c) BFKL (inc. enhanced) 




Figure 4: Sketches of the basic diagram for semi-hard particle production in pp collisions. The 
figure is taken from Ref. [8J. 

On the contrary, the DGLAP-based Monte Carlos generate parton cascades strongly ordered 
in kt- That is, the parton k t increases as we evolve from the input PDF of the proton to the 
matrix element of the hard subprocess, which occurs near the centre of the rapidity interval, 
Fig. |4^b). Since the cross section of the hard subprocess behaves as da/dk^ oc 1/kf, the 
dominant contributions come from near the lower limit k m i n , of the k t integration. In fact, 
in order to describe the high-energy collider data, it is necessary to artificially introduce an 
energy dependent infrared cutoff; k miQ oc s a wittQ a ~ 0.12 [9j. This cutoff is only applied to 

6 This value of a is very close to that obtained from the resummed NLL BFKL prediction that the saturation 
scale satisfies Qg at (x) oc a; -0 45 [SJ. Since x oc 1/y/s and Q sat = fc m in, we expect a = 0.45/4. 
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the hard matrix element, whereas in the evolution of the parton cascade a constant cutoff ko, 
corresponding to the input PDFs, is used. Note that during the DGLAP evolution, the position 
of the partons in 6-space is frozen. Thus such a cascade also forms a hot spot. 

Accounting for the multiple interaction option, that is for contributions containing a few hot 
spots, we include the eikonal multi-Pomeron contributions, both for the DGLAP and BFKL 
based descriptions; that is the presence of a few small-size QCD Pomeron cascades. 

Next, we include the enhanced multi-Pomeron diagrams introducing the absorption of the 
low k t partons. The strength of absorption is driven by the parton density and therefore the 
effect grows with energy, that is with ln(l/x). We thus have an effective infrared cutoff, k sat (x), 
which modifies the kt distribution of the 'BFKL' cascade. The result is shown Fig. Etc), which 
has some similarity to the DGLAP cascade of Fig. 131(b). However, now the cutoff fc sat is not 
a tuning parameter, but is generated dynamically by the enhanced multi-Pomeron diagrams. 
Recall that the same diagrams describe high-mass proton dissociation. That is, the value of 
the multi-Pomeron vertex simultaneously controls the cross sections of high-mass dissociation 
and the effective cutoff fc sat - two phenomena which, at first sight, appear to be quite different. 



4 The Durham model 

Main parameters: 

A = a P (0)-1 , bare Pom. intercept (expect A ~ 0.35) (s-dep.) 

slope a' P (taken zero in the 201 1 KMR analysis) 
d controls BFKL diffusion in log k t 

X strength of triple (and multi) Pom. couplings (triple-Regge data) 

y specifies diff ve estates (determined by low M diff ve dissoc n ) 
N absolute value of initial gluon density 

Figure 5: Some of the main parameters of the KMR model [10]; and how they are constrained. 

How may the partonic model of the Pomeron be implemented in practice? To achieve this 
we note that the absorption of low k t partons is driven by the opacity, Q, which depends both 
on kt and y = m(l/x). The opacity, Qik(y, h, b), is obtained [10] by solving the corresponding 
BFKL-type evolution equations in y with a simplified form of the kernel, but which incorporates 
the main features of BFKL: diffusion in Inkf and A = ojp are (0) — 1 ~ 0.35. (A two-channel 
eikonal is used, i,k — 1,2.) The inclusion of the k t dependence is crucial for the transition 
from the hard to the soft domain. The absorptive factors in the equation embody the result 
that there is less screening for larger kt- The model [10] has only a small number of physically 
motivated parameters, see Fig. [5j whose values are tuned to reproduce the available high 
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energy pp and pp data for cr to t, da e \/dt } a^ M , cr g {f jdtdM 2 etc. Given Q ik (y,k t} b) we can, 
in principle, predict all soft and semi-hard inclusive phenomena, such as the survival factors of 
rapidity gaps, the PDFs and diffractive PDFs at low x and low scales, etc. The predictions for 
some of the cross sections are given in the left-half of Table [TJ as their values will be relevant 
for the discussion in Section [5] It is important to note that hadronization can be incorporated 
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Table 1: Some results of the complete KMR model [TU] prior to the LHC data (left-hand Table), 
and results obtained from a simpler approach, described in Section [5] based on a 3-channel eikonal 
description [llj of all elastic (and quasi-elastic) pp and pp data, including the TOTEM LHC data 
(right-half of the Table). a tot , cr e i and cr^™ are the total, elastic and low-mass single and double 
dissociation cross sections (in mb) respectively. The cross section a SD is the sum of the dissociations 
of both the 'beam' and 'target' protons. B c \ is the mean elastic slope (in GeV -2 ), da e \/dt = e Bclt , 
in the region \t\ < 0.2 GeV 2 . The collider energies are given in TeV. The former (latter) analysis 
fit to the CERN-ISR observations that aj^, A/ =2(3) mb at y/s = 53 GeV, with low mass defined to 
be M < 2.5(3) GeV. 

in this partonic description of the Pomeron, via Monte Carlo generators, which now would have 
the advantage of an effective dynamical cutoff fc sat to suppress low k t parton emissions. 

In summary, some of the main features of the model are: 

(i) values of the high energy pp total cross section which are suppressed by absorptive 
corrections. Increasingly large values of 0"gJf hM are found due to the increasing phase space 
with collider energy. 

(ii) multi-Pomeron contributions arising from eikonal diagrams, that is the presence of a 
few small-size QCD Pomeron cascades (hot spots). This can be tested by measuring Bose- 
Einstein correlations, see Fig. |6j Specifically, identical pion correlations measure the size of 
their emission region. 

(iii) multi-Pomeron contributions arising from enhanced diagrams, which lead to the ab- 
sorption of low k t partons and automatically introduce an effective cutoff k SSut which increases 
with energy. Due to the cutoff, k t > k sat , the main inelastic process is minijet production. The 
dominance of minijets can be tested by observing the two-particle correlations of secondaries 
at the LHC 0. 
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Figure 6: A sketch of the size (r) of the source of identical pions as a function of the multiplicity of 
charged particles in high-energy pp interactions [12]. The continuous and dashed curves correspond 
to, say, yfs = 7 and 0.9 TeV respectively. At low N ch we expect (r) to be independent of collider 
energy, 



r s, while for the plateau we expect (r) ~ R pp oc y B c \(s) to increase very slowly with 
energy. Very high multiplicities are expected to arise from high-.&r events originating from a single 
ladder. 



5 Implications of latest LHC 'soft' data 

As a postscript to the Durham approach, we briefly discuss the implications of some recent 
LHC data on 'soft' diffract ive processes which became available after the LC11 meeting. 

First, we look at the implications of the recent TOTEM measurements, at 7 TeV, of da e \/dt 
down to -t = 0.02 GeV~ 2 [TJj. From these data, TOTEM find 

Otot = 98.3 mb, a c \ = 24.8 mb, <7 inc i = 73.5 mb. (3) 

In the discussion below, we will ignore the (important) experimental errors, just to get some 
ideas of the trends of the data. The KMR model [TU] predicts lower values of 88, 22 and 66 
mb respectively, see the left-half of Table [T} The model was tuned to describe collider data 
for 0"tot- At the Tevatron energy, where the CDF [14] and E710 [15] measurements disagree by 
some 10%, we were much closer to the lower E710 value. 
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To investigate this further, we performed a simpler study than that in [10J. The idea was 
to see if we can describe all the elastic pp and pp collider data in terms of a 3-channel eikonal 
model with only one Pomeron, with parameters that are naturally linked to the perturbative 
QCD (BFKL) framework, as discussed in the previous sections of this paper. However, for the 
simpler study [11], we used an effective Pomeron, rather than the bare QCD Pomeron with 
intercept Ab are = <^p(0) — 1 = 0.32 of [TO]. With an economical parametrization of the three 
(Good- Walker) diffractive eigenstates, we are, indeed, able to obtain a good description of all 
these data for \t\ < 0.3 GeV 2 with A eff = 0.14, see Fig. 7 Since our eikonal model was devised 
to fit the data it is not surprising to have agreement with the TOTEM cross sections of (J3]). 
We call the Pomeron 'effective' since, although we accounted for eikonal rescattering of the 
incoming partons, unlike [TO], we did not explicitly consider enhanced rescattering involving 
intermediate partons. The latter are included implicitly since their main effect is to renormalize 
the bare Pomeron trajectory. 

However, some observations can be made. First, the eikonal description is close to the CDF 
total cross section value, and yet the description of the E710 da e \/dt measurements is good. 
Secondly, the insertion of the pion loopQ into the Pomeron trajectory would decrease the total 
cross value in Table [l] by 1 mb, to 96.4 mb [TT]. The conclusion is that there is quite a bit 
of uncertainty in the extrapolation of the da e \/dt data to the optical point, in addition to the 
3-4% normalization uncertainty. Future precise elastic measurement^] even closer to t = will 
help reduce the uncertainty in the value of cr to t- If the values of cr to t and da e \/dt at the LHC are 
confirmed to be significantly higher than those obtained in [10], then this full analysis should 
be repeated with these data included. It will result in a somewhat larger value of Ab are . 

Let us now compare the results shown in Table [T] with the inelastic cross section obtained 
by CMS, ATLAS and ALICE at 7 TeV. The measured value is defined as the cross section 
with at least two particles in some central (but far from complete) rapidity, r), interval. For 
instance, ATLAS find c>"i nc i = 60.3 mb for the cross section of processes with M > 15.7 GeV, 
that is £ = M 2 /s > 5 x 10~ 6 [21J. After a model dependent extrapolation to cover the entire 
rapidity interval they obtain <r inel = 69.4 mb. CMS find a very similar result, namely 68.0 mb 
[22J. ALICE also get a similar result [23] • These estimates are about 5 mb lower than the recent 
TOTEM value of 73.5 mb of Q. The difference may be attributed to the extrapolated values 
being 5 mb deficient for low-mass diffraction. (The extrapolation in the high-mass interval is 
confirmed by the ATLAS measurement da/dArj ~ da/dhiM 2 ~ 1 mb per unit of rapidity |26j.) 
More specifically, if we define low mass to be M < 3 GeV, then, noting that the unmeasured 
interval from M = 15.7 to M = 3 GeV gives AlnM 2 = 3.3, it follows that the ATLAS, CMS 
results imply o\^ M ~ 64 mb. Then using the TOTEM result we find that low-mass diffractive 

7 Recall that the pion loop insertion modifies the Pomeron trajectory at very small t [IS]. Indeed the presence 
of the 2ir singularity at t = 4mJ leads to some curvature in the t behaviour of dcr e \Jdt. That is, to some variation 
of the local elastic slope B e \(t). Including the pion loop gives an equally good description of the elastic data. 

8 Note also that the simultaneous measurement of bremsstrahlung photons, accompanying elastic proton- 
proton scattering in CMS will, with the help of the Zero Degree Calorimeter, allow an independent determination 
of (Tei/ (B e i); see [20] for details. 
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Figure 7: The description of the data for the differential cross sections forpp (or pp) elastic scattering 
in the energy range 62.5 to 7000 GeV [13, [E3 HU [15] [IT] using a 3-channel eikonal model. The 
Tevatron data with open and closed circles are those of the CDF and E710 collaborations respectively. 
Only very selected TOTEM points are shown, which have been read off their published plot. The 
excellent agreement of our model with the data for small \t\ is achieved with a very economical 
parametrization of the diffractive eigenstates. It is straightforward to describe the elastic data in the 
region of the diffractive LHC dip [TB], but at the expense of a more complicated parametrization of 
the form factors of the three eigenstates. 
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dissociation is expected to have a rather large cross section 




73.5 - 64 



9.5 mb. 



(4) 



Note, however, that the low-mass diffractive dissociation given in Table [TJ 



.SD+DD 
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7.3 + 0.9 



8.2 mb, 



(5) 



in satisfactory agreement with Q. 

Here, it is worth noting that, as seen in [24], by triggering inelastic events with the Tl 
and T2 forward tracking telescopes [25], TOTEM can obtain an independent estimate of the 
total inelastic cross section, analogous to that of the ATLAS measurement [2T] of <r inel , which 
used an information from the MBTS scintillation counters. In particular, note that since 
the T2 telescope extends to larger values of pseudorapidity (5.3 < \r)\ < 6.5), this allows a 
wider coverage of high-mass diffraction (down to M ~ 4 GeV). This, in turn, may allow the 
extrapolation to the total value of er inc i to be performed with reduced uncertainty. 

Another valuable set of soft diffractive measurements have been made by the ATLAS col- 
laboration. They measure da/dArj versus Arj for events with large rapdity gaps [26]. For 
Arj > 5, fluctuations in hadronization are greatly suppressed [27], and we cleanly probe high- 
mass diffractive dissociation. In Ref. [H] these data are shown to be well described by a 
triple-Pomeron approach, provided the sizeable absorptive or rescattering corrections are taken 
into account. These corrections are computed in a parameter-free way using the 3-channel 
eikonal model discussed above. 
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